Physiologia Plantarum

Recent studies have demonstrated the presence of kynurenine (KYN) and kynurenic acid (KYNA) in several plant species, but the metabolic function of these metabolites remains undefined. We hypothesized that KYN and KYNA are metabolites of auxin and play a role in plant morphogenesis. To test our hypothesis, we developed a plant tissue-culture-based bioassay using Hypericum perforatum (St. Johns wort; SJW), a model system for auxin and indoleamine metabolism and pharmacological inhibitors (PF-04859989, RO-61-8048, and KMO inhibitor II, JM6) of human kynurenine pathways enzymes. SJW is an interesting model system because explants root in the absence of plant growth regulators but supplementation of the culture media with 10 M IAA induces a callus response without de novo root organogenesis. Supplementation of the culture media with 10 M KYN increased root number and internodal length relative to basal media. We used a previously validated high-resolution mass spectrometry analytical method to quantify KYN, KYNA, and 3-hydroxyanthranilic acid (3-HAA). KYN, KYNA and 3-HAA were quantified in roots and shoots of SJW grown on basal media. Supplementation of the culture media with 10 M KYN increased the concentration of KYN, KYNA and 3-HAA in roots and shoots. Treatment with 10 M IAA increased KYN and 3-HAA concentration in shoots. Three pharmaceutical candidates that are kynurenine pathway inhibitors in humans were taken up into the tissues from the culture media and increased KYN content as compared to basal control. Together, these data propose a role for KYN in IAA metabolism, shoot and root organogenesis. HighlightsO_LIKynurenine metabolites are detected and accumulate in H. perforatum tissue culture C_LIO_LIIAA redirects metabolism towards accumulation of KYN and 3-HAA in shoots C_LIO_LIExogenous KYN promotes KYNA accumulation C_LIO_LIPharmacological inhibition alters kynurenine pathway metabolite profiles in a tissue-specific manner C_LIO_LIKynurenine and IAA differentially regulate root development C_LI

Soil salinization threatens global agriculture reducing yields, yet the metabolic signals controlling salt-sensitive root plasticity in alfalfa remain unclear. We hypothesize that salinity transiently uncouples the sucrose-trehalose-6-P (Tre6P)- Sucrose non-fermenting kinase 1 (SnRK1) nexus, aligning with a biphasic root metabolic response and altered root architecture. Alfalfa seedlings were grown in a hydroponic system and exposed to 200 mM NaCl, with root samples collected from 1 h to 7 d. While primary root growth and biomass remained unchanged, lateral root development was enhanced under salinity. Early response (1 h-1 d) was characterized by reduced carbon metabolites, low Tre6P, increased malondialdehyde, and SnRK1 activation, with a decline in glycolytic and TCA intermediates. During this phase, sucrose was negatively correlated with both Tre6P and SnRK1. Late response (3-7 d) showed a SnRK1 reactivation, Tre6P recovery, and osmoprotectant accumulation, including increased antioxidant capacity (+75% at 3dpt), proline (+178%), and sucrose (+18%) and starch depletion (-57%) at 7dpt respect to control. These metabolic changes coincided with the enhanced lateral root emergence. These findings indicate a two-phase response: early metabolic downscaling with transient Suc-Tre6P-SnRK1 disruption, followed by recovery with Tre6P restoration, SnRK1 reactivation, osmoprotection, and sustained root plasticity under salinity. HighlightSalinity triggers a temporary metabolic shift in alfalfa roots: plants first conserve energy, then adapt to stress, maintaining lateral root growth and flexible root architecture.

The production of reactive oxygen species (ROS) in response to cadmium (Cd) has been extensively studied, demonstrating that they play a key role in the plants response to this heavy metal. While the role of enzymes like RBOHs has been thoroughly studied, the function of other ROS-producing enzymes, such as peroxisomal glycolate oxidase (GOX), remains largely overlooked. Peroxisomal GOX is a core metabolic enzyme of the photorespiratory pathway occurring in chloroplasts, mitochondria and peroxisomes. Using Arabidopsis (Arabidopsis thaliana) mutants lacking the main peroxisomal GOX genes, GOX1 (gox1-1) and GOX2 (gox2-1) we explored their function in plant response to Cd. Although photosynthetic capacity appears to be affected to the same extent in both mutants under control and Cd stress conditions, GOX2 seems to play a greater role in ROS production in response to the metal. Transcriptomic analyses on WT and gox2-1 pointed to the mitochondrial electron transport chain (mETC) as a target of Cd stress. We further investigated the individual GOX1 and GOX2 functions in mETC regulation and redox state. Although oxidative ratio of mitochondria was higher in both mutants, it was more pronounced in the absence of GOX1. Furthermore, the mETC is affected in both mutants but the regulation of its components differs in each mutant. These results point out the different functions of the two photorespiratory GOX isoforms in Arabidopsis, leading to a better understanding of the photorespiratory pathway.

RNA-seq datasets from medicinal yews are crucial for studying paclitaxel biosynthesis. However, cross-study data analyses are hindered by pronounced batch effects. Here, we compiled 45 RNA-seq samples from three studies across four tissues (bark, leaf, root, stem) and assessed 35 preprocessing pipelines combining six normalization strategies with five batch-effect correction approaches. Unsupervised clustering (HCA, k-means, Grade-of-Membership), evaluated using Jaccard and Adjusted Rand indices, revealed significant variability in batch effect removal. Supervised classification of tissue and project labels (Random Forest and linear/radial SVM) demonstrated improved accuracy in tissue type prediction, highlighting the effectiveness of correction methods. The processed data facilitated the identification of 189 putative ABC transporters across samples, six of which showing a strong correlation to the gene encoding 10-deacetylbaccatin-III-10{beta}-O-acetyltransferase, a key biosynthetic enzyme in the taxol pathway. High expression levels in leaf and bark further support their role in taxane intermediates trafficking in taxol biosynthesis. Structural analysis and molecular docking further supported the selection of these candidates, and the agreement between transcriptomic ranking and docking-based prioritization suggests that these transporters may participate in taxane intermediate recognition, trafficking, or export. These findings demonstrate the importance of normalization and batch effect correction in RNA-seq analysis to advance gene discovery in Taxus species and, more broadly, in plant research. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=152 SRC="FIGDIR/small/723993v1_ufig1.gif" ALT="Figure 1"> View larger version (54K): org.highwire.dtl.DTLVardef@1469162org.highwire.dtl.DTLVardef@1f2c4deorg.highwire.dtl.DTLVardef@15ad821org.highwire.dtl.DTLVardef@123676d_HPS_FORMAT_FIGEXP M_FIG C_FIG

Bread wheat (Triticum aestivum) is a staple food crucial for global caloric intake and food security. The current climate emergency demands the development of sustainable agricultural practices, particularly in the context of drought-induced yield reductions in bread wheat. Microalgae-based biostimulants have emerged as promising tools to enhance crop tolerance to drought stress while concurrently mitigating atmospheric CO2 accumulation. This study characterizes the transcriptomic responses to the foliar application of the microalgae-based biostimulant LRMTM in drought-stressed and fully irrigated wheat plants unveiling its mode of action. Drought stress at the tillering stage significantly altered gene expression activating key pathways related to phosphate starvation response (PSR), inositol phosphate signaling, and tocopherol biosynthesis. The application of the microalgae-based biostimulant LRMTM in drought-stressed plants further enhanced the expression of drought-responsive genes, particularly those involved in PSR and carbon fixation. Specific responses to LRMTM treatment in drought-stressed plants were also found related to abscisic acid (ABA) signaling activating genes involved in stomata closure, which plays a critical role in drought tolerance. In fully irrigated plants, LRMTM treatment was also beneficial modulating circadian rhythms, shade avoidance and attenuating stress responses. Phenotypic analysis showed that LRMTM-treated plants exhibited enhanced drought tolerance, increased height and spike length even under fully irrigated conditions. These results indicate that the microalgae-based biostimulant LRMTM not only enhances wheat response to drought but also promotes growth and productivity in both stressed and non-stressed conditions which could contribute to the development of sustainable agriculture in the face of the current climate challenges.

Cytokinin (CK) N-glucosides are the most abundant CK metabolites in Arabidopsis and most angiosperms, yet their role in cytokinin activity and response is unclear. Here, we examined metabolomic, transcriptomic, and proteomic profiles of seven CK N-glucoside conjugates in detached Arabidopsis leaves across a 144-hour dark-induced senescence (DIS) timecourse. All tested N-glucosides were found to undergo a slow conversion to their corresponding base forms at position-dependent rates, with N9-glucosides releasing base faster than their corresponding N7-glucosides. Conversion during DIS was strictly isoform-specific and not accompanied by coordinated induction of CK biosynthesis genes, arguing against de novo synthesis as the source of accumulated base. Despite progressive base accumulation, N-glucoside-treated leaves produced substantially fewer Differentially Expressed Genes than direct base application at comparable base concentrations, revealing a disconnect between hormone presence and transcriptional output. Unbiased model comparison identified the base:glucoside ratio as a stronger predictor of CK-Two Component Signaling (TCS) gene expression than absolute base concentration, though modulated by base-type-specific receptor affinities. Early proteomic profiling further revealed a coordinated response shared across N-glucosides but largely absent from base treatments. Together, these findings support that CK N-glucosides as kinetically slow, position-dependent reservoirs whose presence in abundance modulate activation of CK-TCS elicited by bioactive forms. HighlightsPhysiology, metabolomic, transcriptomic, and proteomic findings here support CK N-glucosides as kinetically slow, position-dependent reservoirs whose presence in abundance modulate activation of CK-TCS elicited by bioactive forms.

The release and management of pheasants (Phasianus colchicus) in the UK for recreational shooting exerts a range of effects on the ecosystem into which they are released. We studied possible effect of nutrient deposition on epiphytic tree flora at 20 pheasant release sites distributed through England (18) and Wales (2) during winter and spring 2023/24. Sites were all Ancient Semi-natural Woodlands (ASNWs) and had substantial (600-8000 pheasants) in a single release pen. We measured N-sensitive and N-tolerant indicator bryophyte and lichen species on tree trunks near to the pen and then in plots along a transect 100m, 250m, 500m and 1km+ away from the pen. To achieve a gradient of pheasant use, the transects were located in the opposite direction to the game managed / shooting area. We recorded 1.9 times more coverage of N-tolerant lichens and bryophytes combined on selected tree species at the pen-edge compared to the control plots. The relationship showed a decline from the pen edge to 250m away but then stabilised. We also detected higher levels of coverage of N-sensitive tree flora at 100m and 250 m compared to the penedge plot. These measures were also higher at these mid distances compared to the 500m and 1000m plots. We suggest far plots were nearer wood edges and were affected by ambient inputs of aerial N from farmland and other external sources. The overall interpretation is that concentrations of pheasants in and around release pens for several months from late summer until early winter in ASNWs does affect the balance of N-sensitive and tolerant tree flora up to potentially 250m and this is a consideration when locating release pens in and near to sensitive woods.

Micro-dwarf tomato cultivars are increasingly considered for urban and controlled-environment agriculture due to their compact architecture and suitability for high-density planting. However, optimal canopy management strategies for these cultivars remain poorly defined. In this study, we evaluated the effects of different leaf removal intensities on leaf-level physiological performance, fruit yield, and fruit quality in three micro-dwarf tomato cultivars (Mohammed, Hahms Gelbe Topftomate, and Red Robin) grown under contrasting seasonal light conditions. Plants were subjected to low (15%), moderate (30%), or severe (90%) leaf removal, and leaf-level gas exchange was measured across canopy layers, along with yield and fruit quality assessments. Severe leaf removal (90%) increased carbon assimilation, transpiration, and stomatal conductance in middle and lower canopy leaves by up to approximately twofold compared with control plants, indicating improved light availability at the leaf level. However, these physiological enhancements did not consistently translate into higher yield, reflecting reduced whole-plant source capacity under excessive leaf removal. Low to moderate leaf removal (15-30%) generally increased or maintained yield and fruit number, whereas severe leaf removal reduced yield in Hahms Gelbe and Red Robin, particularly under low seasonal radiation. In contrast, Mohammed exhibited yield increases of up to 220% under low leaf removal and maintained increased yield even under severe leaf removal under high-light conditions. Fruit quality was largely unaffected by leaf removal, except for total soluble solids, which declined by approximately 12% under severe leaf removal across cultivars, consistent with sugar dilution under source limitation. Overall, these results demonstrate that optimal leaf removal in micro-dwarf tomatoes requires balancing improved canopy light distribution with maintenance of sufficient leaf area for carbon assimilation. Leaf removal thresholds are strongly cultivar- and light-dependent, emphasizing the need for cultivar-specific canopy management strategies in compact tomato systems and controlled-environment agriculture.

Chilling stress induces photosystem I (PSI) photoinhibition in chilling-sensitive cucumber, in which insufficient activity of the chloroplast NADH dehydrogenase-like complex (NDH) leads to PSI over-reduction and damage. However, it is not yet clear whether these findings can be generalized to other species or what the molecular mechanism underlying impaired NDH function is. In this study, we first examined whether NDH is essential for PSI protection under chilling stress using an NDH-deficient rice mutant. Compared with wild-type plants, the NDH-deficient mutant exhibited enhanced PSI over-reduction and pronounced PSI photoinhibition under chilling stress. In contrast, rice plants expressing flavodiiron protein (FLV), which functions as an alternative electron acceptor downstream of PSI, did not exhibit PSI photoinhibition under chilling stress, demonstrating that electron sink capacity of NDH is important for PSI protection under chilling stress. Furthermore, analysis of the factors responsible for NDH dysfunction under chilling stress in cucumber revealed that chilling stress destabilizes the PSI-NDH supercomplex, leading to NDH monomerization and a consequent loss of NDH activity. This NDH monomerization is likely attributable to chilling-induced damage to the light-harvesting complex Lhca, which mediates the association between PSI and NDH. Together, these results indicate that NDH is essential for protecting PSI from photoinhibition under chilling stress in both rice and cucumber, and that chilling-induced destabilization of the PSI-NDH supercomplex represents a key molecular mechanism underlying PSI over-reduction and photoinhibition.

O_LISymbiosis between legumes and rhizobia is beneficial on nutrient-poor soils, as it enables the fixation of atmospheric N2. To establish this symbiosis, gene expression in both the host plant and the symbiont has to be regulated. To understand the underlying RNA-mediated regulation of host gene expression, we designed experiments to identify competing endogenous networks involving circular RNA, microRNA, and linear transcripts during symbiosis, using wt and symbiosis-deficient Lotus japonicus mutants with the rhizobium Mesorhizobium loti (M. loti). C_LIO_LICircRNA, miRNA, and linear transcripts were identified from Lotus japonicus wildtype and CCamK mutant (ccamk-13; snf-1) seedlings without inoculation or with M. loti inoculation using deep short-read sequencing with rRNA-depletion and random primers. C_LIO_LIDifferentially expressed miRNAs showed negative correlations to predicted target genes and may regulate symbiotic processes. The symbiosis essential iron-sensor LjnsRING/BRUTUS expresses a circRNA which was upregulated in symbiotic treatments. This circRNA may act as a target mimic and contribute to nodule longevity. CircRNAs are predicted to act predominantly as trans-regulatory molecules with similar frequencies in Arabidopsis thaliania, Oryza sativa, and Lotus japonicus. C_LIO_LIWe identified novel miRNAs, long noncoding RNAs, and circRNAs, and nominated several as potential new regulatory non-coding RNAs that may act as target mimics to stabilize genes and support symbiosis. C_LI SummarySymbiosis between Lotus japonicus and Mesorhizobium loti involves treatment-specific regulation of competing endogenous RNA networks involving circular RNA, miRNA, and linear transcripts.

Crop improvement can enhance food security, but side effects, such as trade-offs between valuable agronomic traits, are common. Likewise, fertilisation helps ensure high yields, but can devalue mutualisms with soil microbes that would otherwise be essential for nutrient acquisition. If the need for nutritional mutualisms is reduced in crops, mutualisms could be disrupted by selection relaxation or allocation trade-offs. Thus, crops could achieve high yields in spite of, or because of, disruption of nutritional mutualisms. Alternatively, the highest-yielding plants might flourish because they maximise nutrient acquisition from both symbionts and the soil. Here, enhanced mutualism could evolve over the course of agricultural crop improvement. To investigate whether high yields in cultivars and wild accessions are negatively or positively genetically correlated with outcomes in the legume-rhizobia mutualism, we measured whether yield and symbiosis traits trade-off or are positively genetically correlated among cultivars and wild accessions. We also tested whether this relationship differs between accessions released before or after 1950. We measured genetic correlations between yield and mutualism traits in 87 domesticated pea (Pisum sativum) accessions in a common garden agricultural field across three years. Seed yield and N2 fixation (%Ndfa) were positively genetically correlated. While N fixation was more strongly predictive of yield in the pre-1950 accessions than the post-1950 accessions, the underlying positive genetic correlation between the traits did not differ between the groups. The positive genetic correlation between yield and N2 fixation indicates that selection to increase yields has maintained or increased the benefits of the rhizobial mutualism in pea. Our findings predict that breeding to increase yield may continue to produce pea cultivars that get a greater proportion of their N from rhizobia, enhancing symbiotic mutualism and reducing the proportion of N supplied by fertilisation.

Background and AimsPlants deploy triterpenoid saponins as chemical defences against herbivores, yet it remains unclear whether insect digestion detoxifies these compounds or generates equally or more active metabolites. Because saponin bioactivity depends strongly on glycosylation patterns, we examined the fate and defensive activity of hederagenin-derived saponins during herbivory. MethodsLarvae of Plutella xylostella were fed leaf discs containing structurally defined hederagenin-derived saponins. Saponin composition in treated leaves and larval frass was analysed by LC- qTOF-ESI-MS/MS. Feeding assays were used to compare the antifeedant activity of mono- and bidesmosidic forms. Key ResultsLarvae selectively metabolized complex hederagenin-derived saponins into simpler forms, with cellobiosides converted into monoglucosides during digestion, resulting in a marked shift in saponin composition between ingested material and frass. Feeding assays showed that monodesmosidic saponins strongly deterrer feeding, whereas bidesmosidic saponins were largely inactive. The loss of activity in bidesmosidic saponins was not explained by differential metabolism, indicating that glycosylation patterns directly determine biological function. ConclusionsInsect herbivores selectively modify saponin structures through deglycosylation, thereby altering their defensive properties. Our findings demonstrate that glycosylation governs both saponin activity and metabolic fate, highlighting insect-driven turnover as a critical component of plant chemical defence during plant-herbivore interactions. Issue SectionOriginal article

Pennycress (Thlaspi arvense L.) is an intermediate winter oilseed crop that has only recently been domesticated for agronomic use. Improving agronomic traits requires sources of genetic variation, and mutagenesis is frequently used to help overcome the limitations of natural populations. We investigate the impact of Ethyl methanesulfonate (EMS) on genetically effective cells (GECs) to characterize the intra-individual genetic variation of EMS mutagenesis in pennycress. We identified that pennycress contains at least 4 GECs which, when treated with EMS, create unique mutations across different branches within the same individual plant. We then propagated the M2 plants for whole genome sequencing, providing extensive characterization of the EMS mutation profile and developing a gene index as a resource for future reverse genetic screenings. Article SummaryPennycress is an emerging winter oil seed crop in the American Midwest. Domestication efforts have advanced rapidly through a combination of genetic techniques. One of the most successful methods has been the use of a mutant gene index, a large collection of pennycress seed where new genetic variation has been created through Ethyl methanesulfonate (EMS). EMS mutations are not uniform however, and a single treated seed can have wide genetic variation within the resulting plant. We investigate the role of genetically effective cells on EMS variation, and present the full EMS population as a resource for further pennycress domestication efforts.

Salinity is a major crop production constraint in dry pea (Pisum sativum L.), making the development of salt-tolerant varieties essential to improve crop productivity and land-use efficiency. The genetic mechanisms of salt tolerance in dry pea is largely unknown, and research on salt-tolerant genes is limited. In this study, we established comprehensive genomic and transcriptomic resources, along with a robust screening protocol, to dissect the genetic basis of salinity tolerance using two germplasm sets: the USDA pea diversity panel, consisting of approximately 200 globally sourced accessions, and a set of 300 modern elite lines from the NDSU Pulse Crops Breeding Program. Genetic variation for the salinity response was assessed based on ten phenotypic traits, with root dry weight, shoot dry weight, and specific root length identified as key indicators based on their heritability. Genome-wide association mapping uncovered significant genomic regions and several candidate genes linked to salt stress, with the strongest association found on chromosome 6. Overlapping QTL signals across traits suggest a shared genetic architecture underlying salinity tolerance. Field-based transcriptomic analysis further identified five putative genes involved in salinity response conserved across multiple crop species. Notably, Psat5g000800, encoding a glycosyl hydrolase gene, was markedly upregulated under salinity stress. These findings highlight the complex, multi-gene regulatory nature of salinity tolerance in dry pea and underscore the importance of functional validation of candidate genes. This study provides key insights and practical tools to support breeding efforts aimed at improving salt tolerance in dry pea.

Apples (Malus x domestica) are popular fruits grown in temperate regions of the world. The various genotypes must meet a specific threshold amount of cold exposure before they are competent to break dormancy, a quantity approximated as "chill hours". Several varieties have been identified that exhibit an ultra-low-chill requirement, or more precisely shallow dormancy, breaking vegetative and floral buds early in spring in response to minimal cold exposure. These ultra-low-chill genotypes originated from the Bahamas ( Dorsett Golden,1960s), Israel ( Anna, 1950s) and Alabama, USA ( Shell of Alabama, 1880s). The separation in time and space implies that each would feature distinct genetic lesions that govern dormancy control, providing discrete mechanisms to incorporate a low-chill trait in variety improvement. However, analysis of microsatellites and ultimately genome sequence indicates that Dorsett Golden and Anna share strong concordance with the Shell of Alabama genotype, as well as other ultra-low-chill varieties. Kinship analysis confirms that all are closely related, despite differences in year and place of origin. All three low-chill genotypes share common mutations in the DORMANCY ASSOCIATED MADS-BOX1(DAM1) gene, a known repressor of vegetative growth during dormancy. Genomic sequence diversity is observed among Shell of Alabama individuals, including differences in DAM1 that match differences in flowering time. The results of this study call into question the pedigrees of the ultra-low-chill apple germplasm and indicate variation in an otherwise narrow genetic base for use in future breeding efforts.

RAP2.6, an AP2/ERF transcription factor (TF), regulates plant stress responses; however, its role in floral transition remains unexplored. Here, we evaluated RAP2.6s role in flowering and the associated transcriptional changes in Arabidopsis thaliana under long-day conditions. RAP2.6-overexpressing line showed early flowering with fewer rosette leaves, whereas rap2.6-1 mutant flowered later, had more rosette leaves, and higher expression of the floral repressor FLOWERING LOCUS C (FLC). Early flowering in the overexpressing line was accompanied by transcriptional activation of the floral integrators GIGANTEA (GI), FLOWERING LOCUS T (FT), and COSTANS (CO), potentially through RAP2.6 interaction with GCC/DRE cis-regulatory elements. RAP2.6-mediated floral transition depended on nitric oxide (NO), with flowering time largely varying based on NO bioactivity. RAP2.6 was found to be a downstream regulator of Arabidopsis S-NITROSOGLUTATHIONE REDUCTASE 1 (GSNOR1) in controlling S-nitrosothiol (SNO) levels, flowering time, and silique formation. The NITRIC OXIDE-ASSOCIATED 1 (NOA1)-dependent reduction in NO levels abolished early flowering in 35S::RAP2.6 plants without affecting silique formation. Furthermore, enhanced cytokinin sensitivity and upregulation of cytokinin biosynthetic genes suggest cytokinin involvement in RAP2.6-mediated flowering. Together, these findings highlight the crucial role of RAP2.6 in regulating flowering time by integrating redox and hormonal signaling to coordinate reproductive development in A. thaliana.

The performance of a single cross is determined by the average additive effects of the parents, as well as the interactions between them. These quantities can be estimated using an appropriate genetic design, allowing for the estimation of general (GCA) and specific (SCA) combining abilities. The prediction of GCA for new parents and the total genetic value of unrealized crosses can be made when genome-wide marker information is available. Several studies in crops such as maize and rice have demonstrated the potential of genomic-assisted prediction of single-cross performance in economically important crops. However, no study to date has explored its relevance in perennial ryegrass, an obligate allogamous species that is bred in genetically heterogeneous families. In this study, we aimed to estimate genetic parameters and assess the ability of genomic models to predict the performance of F2 families in terms of dry matter yield and nutritive quality traits. We used data from a large partial diallel involving 104 parents from two distinct subpopulations, as inferred by admixture analysis. F2 families were evaluated in multiple environments and under two nitrogen availability conditions. Genotyping-by-sequencing of the parent plants produced 42,145 variants after quality control, which were used to estimate genomic relationships based on identity-by-state. Variance component estimation revealed limited GCA and SCA interactions with the environment, and particularly with nitrogen management. The predictive abilities of two parental models exceeded 0.60 and often surpassed 0.70 for most traits. However, incorporating non-additive effects into the model did not improve predictive ability. We leveraged the genetic diversity among parents to map genomic regions associated with all recorded traits. Genome-wide association studies (GWAS) by genomic best linear unbiased prediction (GBLUP) identified six quantitative trait loci (QTL) regions, with 45 candidate genes within the linkage disequilibrium range, estimated at approximately 92 kb. Our results demonstrate that genomic prediction of single crosses can be performed with high accuracy, especially when both parents are also progenitors of families in the training set.

(1) BackgroundSaffron crocus (Crocus sativus) is the source of saffron, the most expensive spice in the world. It evolved about 3000 years ago as a sterile triploid clone in Greece. Since then, saffron has spread across the globe, where regionally distinct practices of saffron cultivation have developed. Despite differences in morpho-physiological traits, genetic variability is low, if present at all. Here, we aim to resolve chromosomal and sequence-associated variability across saffron crocus cultivars from the crops main cultivation areas in Africa, Asia and Europe. (2) MethodsWe used genome-wide DNA polymorphisms obtained through genotyping-by-sequencing (GBS) of 33 saffron and 14 closely related Crocus accessions, which we place into a phylogenetic context. For karyotyping, we compare nine saffron accessions by multi-color fluorescent in situ hybridisation (FISH) with repetitive DNA probes. (3) Key resultsPhylogenetic analyses confirmed the single origin and clonal nature of all saffron accessions. We detected slight DNA differences among saffron crocus genotypes, which were minor compared with those in wild C. cartwrightianus populations. Still, the Iranian saffron accessions form a genetically very narrow group that differs from the other proveniences in population genetic analyses. However, chromosomes of some saffron accessions display variable FISH signals, likely resulting from gains and losses of tandemly repeated DNA. (4) Main conclusionsBased on the high genetic identity and small karyotypic differences, we confirm the clonal origin of the saffron accessions. Nevertheless, as we detected small and regional chromosomal variability, we conclude that at least four somaclonal saffron lineages emerged after saffrons origin. Societal Impact StatementFor millennia, many cultures developed cultivation practices and regional crop varieties. A notable case is saffron, the worlds most expensive spice that is harvested from stigmas of saffron crocus. This flower crop arose 3000 years ago in a singular genome triplication event and since then spread clonally across the globe. By identifying genetic and chromosomal variability in clonal saffron accessions, we highlight regional diversity, support the preservation of traditional knowledge, and underscore the risk of relying on only one clonal lineage. This informs strategies for saffron cultivation, linking cultural heritage with modern genomics to address biodiversity, evolution, and food security.

We developed plant (Parallel Annotation of Transcriptomes), a de novo method that can potentially compare RNA-seq data of any two species without a reference genome. plant is conceptually similar to chromatography. In the same way a complex mixture is filtered to isolate its individual components, we applied a computational method to identify, annotate, and quantify components across transcriptomes. The comparison points are universal protein domain annotations rather than species-specific genes, as would be the case for a differential gene expression analysis. We looked at several Selaginella species via the 1000 Plant transcriptomes initiative (1KP) where RNA-seq data for various plant species have been made publicly available. The raw reads were assembled via Trinity. The assembled transcripts were then searched against the Pfam protein domain database via InterProScan. The assembled transcripts were also quantified via kallisto. By merging these two aspects, we were able to see how often a particular protein domain - a predicted protein structure - is expressed. These quantified annotations of protein domains are comparable across species, assuming a relatively short evolutionary distance. We were also able to identify the presence of species-specific protein domains and trace each annotation back to the gene. A bubble plot was created to visualize the distributions of Pfam annotations across species as well as GO terms.

Iron (Fe) is an essential micronutrient for plant growth, and the hormone auxin is a key regulator of developmental processes, including root gravitropism. Here, we investigated the molecular mechanisms underlying the crosstalk between iron nutrition and auxin-mediated root growth in Arabidopsis thaliana. Phenotypic analysis revealed that iron deficiency strongly shaped root system architecture and root gravitropism, and these phenotypes were exacerbated in the iron uptake mutant irt1-1. Genetic analysis revealed that iron deficiency did not aggravate the gravitropic defect of the pin2 mutant, eir1-4, suggesting that iron availability modulates root gravitropism through a PIN2-dependent pathway. Further transcriptomic analysis confirmed that iron deficiency significantly altered the expression of numerous genes related to the auxin pathway, providing molecular evidence for the observed physiological connection. Collectively, this study revealed that iron availability regulates root gravitropic growth by modulating PIN-mediated auxin transport and distribution, providing insights into how plants integrate nutritional cues with developmental programs. Graphical abstract A brief descriptionIron modulates auxin transport and root tip distribution by regulating PIN2 protein, thereby mediating root gravitropism in Arabidopsis. Public summaryO_LIIron nutrition specifically regulates root gravitropism and architecture in Arabidopsis. C_LIO_LIIron deficiency disrupts local auxin homeostasis in root tips and impairs asymmetric distribution in response to gravity. C_LIO_LIIron deficiency stress significantly reduces the abundance of PIN2 protein in root tip cells and disrupts its polar localization pattern on the plasma membrane, thereby precisely modulating polar auxin transport by interfering with the vesicle trafficking and recycling efficiency of PIN2. C_LIO_LIRNA-seq results showed that iron deficiency induced differential expression of multiple auxin-related genes, indicating that iron nutrition affects root development through the auxin pathway. C_LI